Dielectric, piezoelectric, and elastic properties of BaTiOs/SrTiOs ferroelectric 

superlattices from first principles 
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The effect of epitaxial strain on the phonon spectra, crystal structure, spontaneous polarization, 
dielectric, piezoelectric, and elastic properties of (OOl)-oriented ferroelectric (BaTi03) m /(SrTi03) ?l 
superlattices (m — n = 1-4) was studied using the first-principles density-functional theory. The 
ground state of free-standing superlattices is the monoclinic Cm polar phase. Under the in-plane 
biaxial compressive strain, it transforms to tetragonal P4mm polar phase, and under the in-plane 
biaxial tensile strain, it transforms to orthorhombic Amm2 polar phase. When changing the in- 
plane lattice parameter, a softening of several optical and acoustic modes appears at the boundaries 
between the polar phases, and corresponding components of dielectric, piezoelectric, and elastic 
tensors diverge critically. The comparison of the mixing enthalpy of disordered Bao.sSro.sTiOa solid 
solution modeled using two special quasirandom structures SQS-4 with the mixing enthalpy of the 
superlattices reveals a tendency of the BaTi03-SrTi03 system to short-range ordering and shows 
that these superlattices are thermodynamically quite stable. 

PACS numbers: 64.60.-i, 68.65.Cd, 77.84.Dy, 81.05. Zx 



I. INTRODUCTION 

The success in creating of ferroelectric superlattices 
with a layer thickness controlled with an accuracy of 
one monolayer offers new opportunities for design of new 
ferroelectric multifunctional materials with high spon- 
taneous polarization, Curie temperature, dielectric con- 
stant, and large dielectric and optical nonlinearities. Be- 
cause of many problems encountered in the growth and 
experimental studies of ferroelectric superlattices, first- 
principles calculations of their physical properties can be 
used to reveal new promising fields of investigations and 
applications of these materials. 

Earlier studies of thin epitaxial films of ferroelectrics 
with the perovskite structure have shown that their prop- 
erties differ strongly from those of bulk crystals. It 
was established that substrate-induced strain (epitaxial 
strain) has a strong influence on the properties of films. 
Due to strong coupling between strain and polarization, 
this strain changes significantly the phase transition tem- 
perature and can induce unusual polar states in thin 
films-i-i 

To date, the most experimentally studied ferroelec- 
tric superlattice is the BaTi0 3 /SrTi0 3 (BTO/STO) 
one£~— Studies of these superlattices from first princi- 
ple a 24 i 25 ' 33 ~ — have established main factors responsible 
for the formation of their polar structure. The specific 
feature of the superlattice is that the strains induced in 
it by the lattice mismatch between BaTi0 3 and SrTi0 3 
and by the substrate result in concurrency of equilibrium 
polar structures in neighboring layers, so that the po- 
lar structure of the superlattice can be tetragonal, mon- 
oclinic, or orthorhombic, depending on the mechanical 
boundary conditions at the interface with the substrate. 

Although some properties of BTO/STO superlattices 
have been already studied, a number of problems re- 



main unresolved. For instance, first-principles study 
of dielectric properties of these superlattice s 36 ' 37 have 
found only PAmm and Cm polar phases, whereas the 
Ammfl phase, which is characteristic for stretched films 
of BaTi0 3 , 2 SrTi0 3 , 45 and for PbTi0 3 /PbZr0 3 super- 
lattices^ 6 - was not observed. The piezoelectric properties 
were calculated only for PbTi0 3 /PbZr0 3 superlattice^ 
for BTO/STO superlattices these data are absent. Fi- 
nally, the elastic properties of ferroelectric superlattices 
and their behavior at the boundaries between different 
polar phases have not been studied at all. 

In this work, first-principles density-functional cal- 
culations of the phonon spectra, crystal structure, 
spontaneous polarization, dielectric, piezoelectric, and 
elastic properties for polar phases of (OOl)-oriented 
(BTO) m /(STO) n superlattices (SL m/n) with m = n = 
1-4 are performed. The influence of compressive and 
tensile epitaxial strain on the structure and properties of 
polar phases is studied in details for (BTO)i/(STO)i su- 
perlattice. The stability ranges of tetragonal, monoclinic, 
and orthorhombic phases are determined. The critical 
behavior of static dielectric, piezoelectric, and elastic ten- 
sors at the boundaries between different polar phases are 
studied. The ferroelastic type of the phase transitions 
between the polar phases is established. In addition, an 
important question about the thermodynamic stability 
of BTO/STO superlattices is considered. 

The remainder of this paper is organized as follows. In 
Secim we give the details of our calculations. Next, we 
present the results for the ground state (Sec. IIII A[) and 
the polarization (Sec.[nTB]) of (BTO)„/(STO)„ superlat- 
tices. Dielectric, piezoelectric, and elastic properties of 
(BTO )i/ (STO) i superlattice are described in Secs-HHCl 
IIII Dl and IIII El respectively. The thermodynamic stabil- 
ity of the superlattices is analyzed in Sec. IIII Fl The 
obtained results are discussed in Sec. IIVI 
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TABLE I: Parameters used for construction of pseudopotentials.-* Non-relativistic generation scheme was used for Sr, Ti, and 
atoms, and scalar-relativistic generation scheme was used for the Ba atom. All parameters are in Hartree atomic units except 
for the energy V\ oc , which is in Ry. 



Atom 


Configuration 


r s 


r v 


Td 


Qs 


q P 


Id 


^ min 


"max 




Sr 


4s 2 4p 6 4d°5s° 


1.68 


1.74 


1.68 


7.07 


7.07 


7.07 


0.01 


1.52 


1.5 


Ba 


5s 2 5p 6 5d°6s° 


1.85 


1.78 


1.83 


7.07 


7.07 


7.07 


0.01 


1.68 


1.95 


Ti 


3 S 2 3p 6 3d°4 S ° 


1.48 


1.72 


1.84 


7.07 


7.07 


7.07 


0.01 


1.41 


2.65 





2s 2 2p 4 3d° 


1.40 


1.55 


1.40 


7.07 


7.57 


7.07 









II. CALCULATION DETAILS 

The calculations were performed within the first- 
principles density-functional theory (DFT) with pseu- 
dopotentials and a plane-wave basis set as implemented 
in ABINIT software^ The local density approximation 
(LDA) for the exchange-correlation functional^ was 
used. Optimized separable nonlocal pseudopotentials^ 1 - 
were constructed using the OPIUM software^ to improve 
the transferability of pseudopotentials, the local poten- 
tial correction was added according to Ref. [53]. Param- 
eters used for construction of pseudopotentials are given 
in Table |U the results of testing of these pseudopotentials 
and other details of calculations can be found in Ref. HH. 
The plane- wave cut-off energy was 30 Ha (816 eV). The 
integration over the Brillouin zone was performed with 
a 8x8x4 Monkhorst-Pack mesh. The relaxations of 
the atomic positions and the unit cell parameters were 
stopped when the Hellmann-Feynman forces were below 
5 • 10~ 6 Ha/Bohr (0.25 meV/A). 

The lattice parameters calculated using the pseudopo- 
tentials were a = 7.3506 Bohr (3.8898 A) for bulk SrTi0 3 
and a = 7.4923 Bohr (3.9648 A), c = 7.5732 Bohr 
(4.0075 A) for tetragonal BaTiC>3. Slight underestima- 
tion of the lattice parameters (in our case by 0.4-0.7% 
compared to the experimental data) is a known problem 
of LDA calculations. 

The calculations were performed on two structures: 
supercells of 1 x 1 x 2n perovskite unit cells for 
(OOl)-oriented (BTO)„/(STO)„ superlattices (n = 1 
4) and two special quasirandom structures (SQS) for 
Bao.sSro.sTiOa solid solution; the construction of SQSs 
is described in Sec. MI Fl 

Phonon spectra, dielectric, piezoelectric, and elastic 
properties of the superlattices were calculated within the 
DFT perturbation theory. Phonon contribution to the di- 
electric constant was calculated from phonon frequencies 
and oscillator strengths^ The Berry phase method 55 
was used to calculate the spontaneous polarization P s . 

As the layers of the superlattices are epitaxially grown 
on (OOl)-oriented substrate with a cubic structure, the 
calculations were performed for pseudotetragonal unit 
cells in which two in-plane translation vectors have the 
same length ao and all three translation vectors are per- 
pendicular to each other. This means that for monoclinic 
and orthorhombic phases (with Cm and Amm2 space 



groups) small deviations of the angles between the trans- 
lation vectors from 90° (which were typically less than 
0.07°) were neglected. As was checked, this does not 
influence much the results. 



III. RESULTS 

A. Ground state of epitaxially strained 
superlattices 

The lattice mismatch between BaTiC>3 and SrTiC>3 
creates tensile biaxial strain in SrTi03 layers and com- 
pressive biaxial strain in BaTi03 layers of free-standing 
BTO/STO superlattice. If the layers were isolated, these 
strains would result in appearance of the in-plane sponta- 
neous polarization in SrTiOs (Amm2 space group) and in 
increasing of the out-of-plane polarization in BaTi03 lay- 
ers (Cm or PAmm space groups)^— As the polar state 
with strong local variations of polarization is energet- 
ically unfavorable^ the structure of the polar ground 
state of the superlattice requires special consideration. 
Earlier studies of BTO/STO^^ and PbTiOs/PbZrO^ 
superlattices have demonstrated that both the magni- 
tude and orientation of polarization depend also on the 
substrate-induced (epitaxial) strain in superlattices. 

The ground state of BTO/STO superlattice was 
searched as follows. For a set of in-plane lattice param- 
eters ao, which were varied from 7.35 to 7.50 Bohr, we 
first calculated the equilibrium structure of the paraelec- 
tric phase with PA/mmm space group by minimizing the 
Hellmann-Feynman forces. The phonon frequencies at 
the r point were then calculated for these structures. It 
is known that the ground state of any crystal is character- 
ized by positive values of all optical phonon frequencies 
at all points of the Brillouin zone. So, if the structure 
under consideration exhibited unstable phonons (with 
imaginary phonon frequencies), the atomic positions in 
it were slightly distorted according to the eigenvector of 
the most unstable mode, and a new search for the equi- 
librium structure was initiated. The phonon frequencies 
calculation and the search for equilibrium structure were 
repeated until the structure with all positive phonon fre- 
quencies was found. 

It should be noted that the only unstable mode in the 
paraelectric PA/mmm phase of (BTO)i/(STO)i super- 
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FIG. 1: Phonon spectra for (a) the PAmm phase of 
(BTO)i/(STO)i superlattice grown on SrTi0 3 substrate 
(do = 7.3506 Bohr) and (b) the Cm phase of the same, free- 
standing superlattice (ao = 7.4461 Bohr). 



lattice is the ferroelectric one at the T point. The well- 
known structural instability of SrTi03 associated with 
the i?25 phonon mode at the boundary of the Brillouin 
zone disappears in the superlattice: the frequency of the 
corresponding phonon at the M point of the folded Bril- 
louin zone (to which the R point transforms when dou- 
bling the c lattice parameter) is 55 cm -1 for 1/1 free- 
standing superlattice and 61 cm -1 for 1/1 superlattice 
grown on SrTiC>3 substrate (see Fig. [TJ). 

The phonon spectra calculations show that in 1/1 su- 
perlattice grown on SrTiC-3 substrate (compressive epi- 
taxial strain, the in-plane lattice parameter ao is equal 
to that of cubic strontium titanate) the tetragonal po- 
lar phase with PAmm space group is the ground state 
(Fig.Q}. For free-standing superlattice, the PAmm struc- 
ture is unstable and transforms to monoclinic Cm polar 
one. Under tensile epitaxial strain (ao = 7.46 Bohr), the 
orthorhombic Amm2 polar phase is the most stable one 
for 1/1 superlattice. This means that the variation of 
ao (for example, by growing the superlattice on different 
substrates) can be used to control the polar state of the 
superlattice. 

In order to determine accurately the location of the 
boundaries between PAmm and Cm phases and between 
Cm and Amml phases for 1/1 superlattice, the ground 
state was calculated for a set of in-plane lattice parame- 
ters ao, and for each of these structures the phonon fre- 
quencies at the r point were computed. The dependence 
of four lowest phonon frequencies on the ao parameter is 
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FIG. 2: (Color online) Frequencies of four lowest phonon 
modes at the F point for polar phases of (BTO)i/(STO)i su- 
perlattice as a function of the in-plane lattice parameter ao- 
The labels near the points indicate the symmetry of modes. 
Vertical lines indicate the phase boundaries. 



plotted in Fig. [2j It is seen that the frequency of a doubly 
degenerate E mode decreases critically when approach- 
ing the boundary between PAmm and Cm phases from 
the tetragonal phase. After transition to the monoclinic 
phase two non-degenerate A' and A" soft modes appear 
in the phonon spectrum; the first of these modes becomes 
soft again when approaching the boundary between Cm 
and Amml phases. The soft phonon mode in the Amml 
phase has the B\ symmetry. 

Extrapolation of the squared frequencies of soft ferro- 
electric modes (E mode in the PAmm phase, A' mode 
in the Cm phase and B\ mode in the Amm2 phase) as 
a function of ao to zero gives the in-plane lattice pa- 
rameters corresponding to the boundaries between dif- 
ferent polar phases. The PAmm-Cm boundary is at 
ao = 7.4023 Bohr when extrapolating from the tetrag- 
onal phase and at 7.4001 Bohr when extrapolating from 
the monoclinic phase. The Cm-Amm2 boundary is 
at ao = 7.4489 Bohr when extrapolating from the or- 
thorhombic phase and at 7.4483 Bohr when extrapolat- 
ing from the monoclinic phase. Small difference between 
the values obtained from extrapolation from two sides 
of the boundary means that both phase transitions are 
close to the second-order ones. Taking into account that 
a zero in-plane strain in polar superlattice corresponds 
to the in-plane lattice parameter of ao = 7.4462 Bohr, 
we get the values of —0.605% and +0.032% for the misfit 
strains corresponding to the phase boundaries. 

Comparison of the energies of different polar phases 
and calculation of the phonon frequencies at the T point 
for (BTO) 2 /(STO) 2 and (BTO) 3 /(STO) 3 superlattices 
shows that the same PAmm, Cm, and Amm2 phases 
has the lowest energy, respectively, in superlattices grown 
on SrTiC>3, free-standing superlattices, and superlattices 
grown on the substrate with ao = 7.46 Bohr. 
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FIG. 3: The in-plane (P±) and out-of-plane (P z ) components 
of polarization for (BTO)i/(STO)i superlattice as a function 
of the in-plane lattice parameter do. Vertical lines indicate 
the phase boundaries. 



FIG. 4: Eigenvalues of the static dielectric constant tensor 
Eij for (BTO)i/(STO)i superlattice as a function of the in- 
plane lattice parameter ao- Vertical lines indicate the phase 
boundaries. 



B. Spontaneous polarization 

The calculated spontaneous polarization for free- 
standing and substrate-supported superlattices, tetrag- 
onal BaTi03, and disordered Bao.sSrrj.sTiOs solid solu- 
tion modeled using SQS-4 structures (see Sec. IIII Fp are 
given in Table [TT1 

As was established in Sec. IIII Al the tetragonal PAmm 
phase with the polarization vector normal to the layers 
is the ground state for BTO/STO superlattices grown on 
SrTiC>3 substrate. The calculations shows that in these 
superlattices the spontaneous polarization P s increases 
monotonically from 0.277 C/m 2 to 0.307 C/m 2 as the 
layer thickness is increased from n = 1 to 4 unit cells 
(see Table [TTJ) . The obtained P s values agree well with 
the value of 0.28 C/m 2 estimated from the data of Ref.[H 
for the superlattice with equal thickness of BaTiC>3 and 
SrTiC-3 layers; the P s value of 0.259 C/m 2 for tetrag- 
onal BaTiOs agrees well with the value of 0.250 C/m 2 
reported in Ref. HH. As follows from Table UH for all 
superlattices the P s values are larger than those for 
Bao.sSro.sTiOs solid solution; for superlattices grown on 
SrTi03 substrates they are even larger than P s of tetrag- 
onal BaTiOs. These results agree with experiment^ and 
results of previous calculations! 33 ' 42 

For free-standing superlattices, the monoclinic Cm 
phase is the ground state; the components of the polar- 
ization vector in this phase are also given in Table [XT] It is 
seen that the polarization is rotated continuously in the 
(110) plane and the magnitude of polarization increases 
with increasing n. 

For biaxially stretched superlattices, the Amm2 phase 
is the ground state and the polarization vector is ori- 
ented along [110] direction of the reference tetragonal 
PA/mmm structure of the paraelectric phase. 

The in-plane and out-of-plane components of the po- 
larization for 1/1 superlattice are plotted as a function 



of the in-plane lattice parameter ao in Fig. [3] Extrap- 
olation of the and P 2 dependence on ao to zero 
gives the positions of the PAmm-Cm and Cm-Amml 
phase boundaries. Their values, ao = 7.4018 Bohr and 
ao = 7.4492 Bohr, are very close to those obtained in 
Sec. IIII Al from the frequencies of soft modes. 



C. Dielectric properties 

The eigenvalues of the static dielectric constant tensor 
= 1, 2, 3) for (BTO)i/(STO)i superlattice as a 
function of the in-plane lattice parameter ao are shown in 
Fig. |U In the tetragonal phase, the eigenvectors of the ex- 
tensor coincide with crystallographic axes and £\\ = £22- 
So, the dielectric properties of this phase are described 
by two nonzero independent parameters, En and £33. 

In the monoclinic phase, the polarization vector rotates 
monotonically in the (110) plane; all three eigenvectors 
of the Sij tensor are different and do not coincide with 
crystallographic axes of the reference tetragonal struc- 
ture. As all nine components of the tensor in this 
coordinate system are nonzero for the Cm phase, the 
most compact way to describe the properties of this ten- 
sor is to present its eigenvalues. In the Cm phase, the 
direction of the eigenvector corresponding to the smallest 
eigenvalue is close, but do not coincide with the direction 
of the polarization. 

In the orthorhombic phase, the eigenvectors of the e,-j 
tensor are oriented along the [110], [110] and [001] direc- 
tions of the reference tetragonal structure. The direction 
of the eigenvector corresponding to the smallest eigen- 
value coincides with the polarization vector and the di- 
rection of eigenvector corresponding to the largest eigen- 
value is [001]. 

As follows from Fig.UJ at least one of the eigenvalues of 
the Sij tensor diverges critically at the PAmm— Cm and 
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TABLE II: Spontaneous polarization (in C/m 2 ) for BTO/STO superlattices with different thickness of layers, two SQS-4 
structures used for modeling of disordered Bao.5Sro.sTi0 3 solid solution, and tetragonal barium titanate. The in-plane lattice 
parameters oo (in Bohr) for superlattices are also presented. 



Structure 


SL1/1 


SL2/2 


SL3/3 


SL 4/4 


SQS-4a 


SQS-4b 


BaTi0 3 


P s orientation 


[xxz] 


[001]° 


[xxz] 


[001]° 


[xxz] 


[001] " 


[xxz] 


[001]" 


[111] 


[001] 


[001] 


P z 


0.061 


0.277 


0.113 


0.293 


0.146 


0.302 


0.157 


0.307 


0.130 


0.206 


0.259 


P — P 


0.165 





0.159 





0.150 





0.144 





0.130 








a 


7.4461 


7.3506 


7.4432 


7.3506 


7.4403 


7.3506 


7.4391 


7.3506 









"Superlattices grown on SrTi03 substrate. 



TABLE III: Largest piezoelectric moduli for monoclinic phase 
of free-standing (BTO)i/(STO)i superlattice, for tetragonal 
phase of the same superlattice grown on SrTiOs substrate, 
and for tetragonal barium titanate. 



Structure 


SL1/1 


BaTiOs 


P a orientation 


[xxz] [001]" 


[001] 


e 33 , C/m 2 (d 33 , pC/N) 


31.9 (460) 7.1 (49) 


6.3 (42) 


eis, C/m 2 (dis, pC/N) 


-0.09 (-19) 3.2 (31) 


-2.9 (-24) 



"Superlattice grown on SrTi0 3 substrate. 



Cm-Amm2 boundaries as the in-plane lattice parame- 
ter a is changed. When approaching the PAmm-Cm 
boundary from the tetragonal phase, the en = £22 com- 
ponents of this tensor diverge as the polarization vector 
P s || [001] becomes less stable against its rotation in the 
(110) plane. When approaching the Cm-Amm2 bound- 
ary from the orthorhombic phase, the £33 value diverges 
as the polarization vector P s || [110] becomes less stable 
against its rotation in the same plane. 



D. Piezoelectric properties 

Due to high sensitivity of both magnitude and orienta- 
tion of the polarization vector in superlattices to epitax- 
ial strain, we can expect them to be good piezoelectrics. 
It is known that anomalously high piezoelectric mod- 
uli found in some ferroelectrics like PbZri-^Ti^Os near 
the morphotropic phase boundary are due to the ease 
of strain-induced rotation of polarization in the interme- 
diate monoclinic phasei 56 ' 57 A similar situation appears 
in the monoclinic phase of BTO/STO superlattice. To 
our knowledge, the piezoelectric properties of BTO/STO 
superlattices have not been studied so far neither exper- 
imentally, nor theoretically. The only superlattice, for 
which some piezoelectric properties were calculated, is 
the PbTi0 3 /PbZr0 3 1/1 superlattice^ which was used 
to simulate the properties of PbTio.5Zro.5O3 solid solu- 
tion. 

The largest piezoelectric stress moduli (i = 1, 
2, 3; v = 1-6) calculated for the PAmm phase of 
(BTO)i/(STO)i superlattice grown on SrTi0 3 substrate 
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FIG. 5: (Color online) Components of the piezoelectric ten- 
sor eiv in polar phases of (BTO)i/(STO)i superlattice as a 
function of the in-plane lattice parameter ao- Vertical lines 
indicate the phase boundaries. 



and for the Cm phase of the same free-standing super- 
lattice are given in Table Mil It is seen that in tetragonal 
phases of the superlattice and BaTi0 3 the e 33 moduli 
do not differ much. However in the monoclinic phase, 
which is the ground state for free-standing superlattice, 
the 633 value is five times larger. Even stronger effect 
can be seen for the ^33 piezoelectric strain coefficient 
{diu = X)u=i e i[iS[iv), which is a result of an 1.5-fold 
increase in the elastic compliance modulus S33 in the 
monoclinic phase (see Sec. IIII E| . 

The piezoelectric moduli ei u in polar phases of 
(BTO)i/(STO)i superlattice as a function of the in-plane 
lattice parameter are shown in Fig. [3] According to the 
symmetry, in the tetragonal phase the piezoelectric ten- 
sor has three independent and five nonzero components: 
e 3i = e 32, 633, and ei5 = e24. In our superlattice only 
two of them have large values: 633 and ei5. When ap- 
proaching the PAmm-Cm boundary from the tetrago- 
nal phase, the 633 value increases monotonically whereas 
the eis value diverges critically and reaches the value of 
80 C/m 2 (not shown). 

In the orthorhombic phase (in the coordinate system 
of the reference tetragonal structure) the en = e22, 
ei2 = e 2 i, ei 3 = e 23 , e 34 = e 35 , and e 16 = e 26 moduli 
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are nonzero, and the total number of independent pa- 
rameters is 5. Among them the en, ei2, 634, and eig 
moduli have the largest values (see Fig. [5]). The only 
modulus that behaves critically at the boundary between 
Cm and Amm2 phases is the 634 one, with a maximum 
value reaching 192 C/m 2 (not shown). 

The most complex behavior of piezoelectric moduli is 
observed in the monoclinic phase because the change of 
the in-plane lattice parameter ao results in monotonic ro- 
tation of the polarization vector in the (110) plane. In 
the coordinate system of the reference tetragonal struc- 
ture, all 18 components of the e^ tensor are nonzero (the 
total number of independent parameters is 10). As fol- 
lows from Fig. [5j in the monoclinic phase the en = e22, 
eis = e 2 4, ei2 = e 2 i, ei 3 = e 23 , ei 4 = e 25 , and e w = e 26 
moduli diverge critically at the boundary between Cm 
and PAmm phases, and 633, 634 = 635, e3i = 632, and 
e36 moduli diverge critically at the boundary between 
Cm and Amml phases. It should be noted that when 
crossing the boundaries, small additional jumps (about 
~10%) are observed in the 633 modulus (at the Cm- 
PAmm boundary) and in en and ei2 moduli (at the 
Cm-Amm2 boundary). 



E. Elastic properties 

It is known that ferroelectric phase transitions between 
two polar phases are often the improper ferroelastic ones, 
which means that they are accompanied by appearance 
of soft acoustic modes and spontaneous strain, but the 
strain is not the primary order parameter. This occurs 
when the strain tensor and the polar vector transform 
according to the same irreducible representation of the 
high-symmetry phased As such phase transitions occur 
in BTO/STO superlattices when changing the in-plane 
lattice parameter ao , it was interesting to study the influ- 
ence of these phase transitions on the elastic properties 
of superlattices, especially taking into account that these 
properties of superlattices have not been studied so far. 

In the tetragonal PAmm phase, the elastic compliance 
tensor S^u (Mj v = 1-6) has six independent and nine 
nonzero components. In the orthorhombic Amml phase 
(in the coordinate system of the reference tetragonal 
structure) the tensor has 9 independent and 13 nonzero 
components, and in the monoclinic Cm phase it has 
13 independent and 21 nonzero components. 

The components of the elastic compliance tensor S^ v 
for polar phases of (BTO)i/(STO)i superlattice are plot- 
ted as a function of the in-plane lattice parameter ao in 
Fig. [6] It is seen that at the boundary between PAmm 
and Cm phases the components of S^ u tensor exhibit a 
step-like change (5 i3 = 523, S33, S 66 , S w = S 2 e, S 36 , 
S34 = S35, S46 = S56 moduli), a critical divergence 
from the monoclinic side (Su — S22, 5*12, Si 5 = S24, 
Su = S25, S45 moduli), or critical divergences from both 
sides of the boundary (S44 = S55 modulus) . In the mon- 
oclinic phase, the S14, S15, Sie, 534, 536, S45, and £46 
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FIG. 6: (Color online) Components of the elastic compliance 
tensor S^v for polar phases of (BTO)i/(STO)i superlattice 
as a function of the in-plane lattice parameter ao- Vertical 
lines indicate the phase boundaries. 



moduli become nonzero. In the orthorhombic phase, the 
S14, S15, 534, and £46 moduli vanish again whereas the 
other moduli remain nonzero. At the boundary between 
Cm and Amml phases the anomalies in elastic moduli 
are smaller: the Su, S12, and Sqq moduli exhibit step- 
like changes, the 513, Su, Sis, Si 6 , S 33 , S34, S^, and 
£46 moduli exhibit weak divergence from the monoclinic 
side, and the £44 and S45 moduli exhibit weak divergence 
from both sides of the boundary. 

The critical divergence of the S44 modulus at both 
PAmm-Cm and Cm- Amml boundaries indicates that 
the phase transitions induced in the superlattice by the 
increase of the in-plane lattice parameter are indeed the 
ferroelastic ones. 



F. Thermodynamic stability 

Thermodynamic stability of ferroelectric superlattices 
is very important for their possible applications. Ther- 
modynamic stability of BTO/STO superlattice is deter- 
mined by the mixing enthalpy of the superlattice and 
its relationship with the mixing enthalpy of the disor- 
dered Bao.sSro.sTiOs solid solution. The most complex 
part of the first-principles calculation of thermodynamic 
stability is the calculation of the mixing enthalpy for a 
solid solution because its simulation using supercells with 
a large number of randomly distributed atoms makes it 
extremely time-consuming. 

A conceptually new approach to this problem was pro- 
posed by Zunger et al£Q In this approach, a disordered 
solid solution A x Bi- x is modeled using a special quasir- 
andom structure (SQS) — a short-period superstructure, 
whose statistical properties (numbers of Naa, Nbb, and 
Nab atomic pairs in few nearest shells) are as close as 
possible to those of ideal disordered solid solution (at the 
same time the sites of the superstructure are determin- 
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TABLE IV: Translation vectors, superlattice axis, stacking sequence of atomic planes, and correlation functions il2, m for SQS-4 
structures used for modeling of disordered Bao.sSro.sTiOa solid solution. 

Structure Translation vectors Axis and stacking sequence Il2,i 112,2 n2,3 112,4 

SQS-4a [211], [112], [121] [111] AABB -1 

SQS-4b [210], [210], [001] [120] AABB -1/3 1/3 




• A 

O B 



a 2 a } 



FIG. 7: (Color online) The unit cell of SQS-4b superstruc- 
ture contains four primitive cells of the perovskite structure 
and is constructed using ai, a-2 and a3 translation vectors. 
The atoms of one of two types, A or B, occupy the sites of 
one sublattice of the perovskite structure lying on the planes 
(shown by blue), which are perpendicular to the superlat- 
tice axis (shown by red arrow). The stacking sequence of the 
planes along the superlattice axis is AABB. 

istically filled with A and B atoms). This method has 
been widely used to study the electronic structure and 
physical properties of semiconductor solid solutions and 
the ordering phenomena in metal alloys. To study the 
properties of ferroelectric solid solutions this approach 
was used quite rarej 61 ' 62 

The structure of the disordered Bao.sSro.sTiOa solid 
solution was modeled using two special quasirandom 
structures SQS-4 constructed with the gensqs program 
from ATAT toolkit.— One of these structures is sketched 
in Fig. [7] The translation vectors, superlattice axis, 
and stacking sequence of the planes filled with the same 
atoms, Ba or Sr (denoted by A and B), are given in Ta- 
blc lIVI The pair correlation functions Il2,m (m is the shell 
number) , which describe the deviation of statistical prop- 
erties of these SQSs from those of an ideal solid solution, 
are also given in this table. For x = 0.5 U2, m is simply 
(2NAA/N m — 1), where N rn is a number of neighbors in 
the mth shell. As follows from this table, for the SQS- 
4a structure strong deviation from an ideal solid solution 
appears only in the forth shell; for the SQS-4b structure 
deviations appear in the second and fourth shells, but are 
smaller. The mixing enthalpy AH for all studied struc- 
tures X (superlattices with different periods and SQS 
structures) was calculated using the formula 

AH = E tot (X) - [£ tot (BaTi0 3 ) + E tot (SrTi0 3 )]/2 

from the values of the total energy E tot (per five-atom 
formula unit) for free-standing fully relaxed paraelectric 
PmZm and PA/mmm phases. The obtained values of 



TABLE V: The mixing enthalpy (in meV) for five 
(BTO)„/(STO) n superlattices with different periods and 
two SQS-4 structures used for modeling of disordered 
Bao.sSro.sTiOs solid solution. 
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13.4 


16.8 


11.0 



AH for these structures are given in Table [V] 

An unexpected result of our calculation is the fact that 
AH values for two shortest-period superlattices (1/1 and 
2/2) appeared smaller than AH values for both realiza- 
tions of disordered solid solution. This means that a 
tendency to short-range ordering of components exists 
in the BaTi03-SrTiC>3 system. Low values of AH for 
these superlattices (< 9 meV) indicate that short-period 
BTO/STO superlattices are thermodynamically stable at 
300 K. 

The tendency to short-range ordering found in (001)- 
oriented BTO/STO superlattices can be explained by 
a general tendency of the A cations to order in a lay- 
ered manner in perovskites, in contrast to the B cations, 
which prefer a rock-salt orderingjSi One can add that 
a similar effect was observed in our studies of (001)- 
oriented (PbTi03)„/(SrTi03) n superlattices, where neg- 
ative values of AH for n = 1-3 and positive values for 
larger n were observed. 

Our values of the mixing enthalpy for BTO/STO su- 
perlattices are much smaller than AH value obtained in 
Ref. H! (42 meV per formula unit). Analysis of the cal- 
culation technique used in Ref. |65| shows that atomic po- 
sitions in the superstructures were not relaxed and the 
superstructures were assumed to be cubic when calculat- 
ing the energies of different atomic configurations. So, 
the calculated mixing enthalpy in this paper includes a 
large energy of excess strain. 



IV. DISCUSSION 

Our results on the influence of epitaxial strain on the 
ground state of BTO/STO superlattice agree o nly par- 
tially with the results obtained in Refs. mUlilll The 
results coincide in that: 1) the ground state for free- 
standing (BTO)i/(STO)i superlattice is the monoclinic 
Cm phase, 2) under the compressive strain, the superlat- 
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tice undergoes the phase transition from Cm to PAmm 
phase, and 3) the dielectric constant diverges at the 
PAmm-Cm phase boundary. At the same time, in con- 
trast to the results of Refs. l36H37l we succeeded to observe 
the phase transition to the Amm2 orthorhombic phase 
under the tensile strain. (In Refs. I36ll37l only the rota- 
tion of the polarization vector towards the [110] direction 
was observed under the tensile strain.) We consider our 
results to be more reliable because they agree with the 
results obtained for strained BaTiOa 2 and SrTiC>3 45 thin 
films, the results of recent atomistic calculations of the 
strain-temperature phase diagram for (BTO)2/(STO)2 
superlattice and with results obtained for another su- 
perlattice, PbTiOs/PbZrOs. 4 — In all these systems the 
same phase sequence, PAmm-Cm-Amm2, was observed 
as the in-plane lattice parameter was increased. 

The increase in the polarization P s in (BTO)„/(STO)„ 
superlattices grown on SrTi03 substrates with increasing 
the layer thickness n (Table [TTJ) agrees with the results 
of Ref. [42 in which the explanation of this phenomenon 
was proposed. In free-standing superlattices, the P z com- 
ponent of polarization also increased with increasing n, 
but the P x and P y components decreased with increas- 
ing n, in contrast to the results observed for 3/3, 4/4 
and 5/5 superlattices with fixed in-plane lattice parame- 
ter equal to 1.01 times the lattice parameter of SrTiC>3. 34 
We attribute these changes to the decrease of the in-plane 
lattice parameter ao for free-standing superlattices with 
increasing n (see Table [n} . 

Unfortunately, only a few data points presented in 
Ref. [13 for the dielectric constant at the PAmm-Cm 
boundary for (BTO)i/(STO)i superlattice did not en- 
abled us to make detailed comparison between our re- 
sults. However, the comparison of our dielectric data 
with those calculated for (PbTi03)i/(PbZr03)i super- 
lattice 4 ^ shows that in the monoclinic phase all eigen- 
values of the £jj tensor for BTO/STO superlattice are 
higher, and so this superlattice may be more promising 
for different applications. 

To obtain large piezoelectric moduli necessary for 
technological applications, the epitaxial strain in the 
BTO/STO superlattice should be tuned to a value at 
which ei v and properties of the superlattice diverge. 
As was shown in Sees. IIII Dl and MI E( this appears at 
the phase boundaries. The calculations shows that in the 
vicinity of the PAmm-Cm boundary the d\i piezoelectric 
coefficient reaches a maximum value of 2300 pC/N in the 
monoclinic phase and the dig coefficient reaches a value 
of 6200 pC/N in the tetragonal phase. In the vicinity 
of the Cm-Amm2 boundary, the ^33 piezoelectric coeffi- 
cient reaches a value of 920 pC /N in the monoclinic phase 
and the ^34 coefficient reaches a value of 10500 pC/N 
in the orthorhombic phase. For comparison, the maxi- 
mum piezoelectric coefficient obtained experimentally on 
single crystals of the Pb(Zn 1 / 3 Nb 2 /3)03-PbTi03 system 
was 2500 pC/N. 66 Anomalous increase in the calculated 
dis and coefficients up to ^8500 pC/N was predicted 
for hydrostatically stressed PbTiOs in the vicinity of the 
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FIG. 8: Frequencies of acoustic modes in the vicinity of the 
r point in polar phases of (BTO)i/(STO)i superlattice as a 
function of the in-plane lattice parameter ao. Vertical lines 
indicate the phase boundaries. 



PAmm-Cm phase boundary^ 

Consider now the elastic properties of BTO/STO su- 
perlattice and the results indicating the appearance of 
improper ferroelastic phase transitions. According to 
Ref. |59L in crystals with PAmm space group the phase 
transition 4mm — > m should be of the ferroelastic type. 
The spontaneous strain at this phase transition is char- 
acterized by one or both nonzero and M5 components 
of the strain tensor, and a soft transverse acoustic (TA) 
mode with a wave vector parallel to the polar axis should 
appear in the phonon spectrum in the vicinity of the 
r point. Direct calculations of frequencies of acoustic 
modes in the tetragonal phase at the point with a re- 
duced wave vector q = (0,0,0.05) confirmed this (see 
Fig. [8]): when approaching the PAmm-Cm boundary, 
the frequency of a doubly degenerate TA phonon with 
the A5 symmetry decreased by 5 times. The softening 
of this mode is directly related to the divergence of the 
S44 = £55 components of the elastic compliance tensor. 

At the boundary between Cm and Amm2 phases, 
there should be another ferroelastic phase transition. Ac- 
cording to Ref. [59l . the transition mm2 — > m is accompa- 
nied by spontaneous strain with one nonzero of three 1*4, 
U5, uq components. In our coordinate system with un- 
usual for orthorhombic phase orientation of the polar axis 
(along the [110] direction) there should be a softening of 
TA phonon with a wave vector oriented along this axis. 
This was confirmed by direct calculations of frequencies 
of acoustic modes at the point in the Brillouin zone with a 
reduced wave vector q = (0.035,0.035,0) (see Fig. [8]). In 
contrast to the tetragonal phase in which the soft acous- 
tic mode is doubly degenerate, in the orthorhombic phase 
the only phonon mode that softens at the phase bound- 
ary is the TA mode polarized along the [001] axis. In our 
opinion, this difference is the reason why the anomalies 
in the elastic properties at the ferroelastic Cm-^4mm2 
phase transition are much weaker than those at the fer- 
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roelastic Pimm-Cm phase transition. Unusual orienta- 
tion of the polar axis in our coordinate system results 
in coupling of some components of the elastic compli- 
ance tensor: for example, the S44 and S45 moduli, which 
diverge in the orthorhombic phase, satisfy the relation 
S44 — 545 ks const in this phase. 

Negative value of the acoustic phonon frequency in the 
Cm phase, which is seen in a narrow wave vector region 
in Fig. [TJ is an artifact of calculations. Computation 
of the phonon dispersion curves in the vicinity of the T 
point revealed three acoustic branches to(q), whose fre- 
quencies increased monotonically with increasing q, but 
gave negative values of u>(0) (about 8 cm -1 ) in the limit 
q — > because of numerical errors. After application of 
the acoustic sum rule, w(0) restored its zero value, but 
the derivative duj/dq near q = for the softest mode 
became negative. So, there is no contradiction between 
the phonon spectra and the positive definiteness of the 
elastic moduli matrix dj calculated in Sec. MI El 

V. CONCLUSIONS 

In this work, the properties of (OOl)-oriented 
(BaTi03) m /(SrTi03)„ superlattices with m = n = 
1-4 were calculated using the first-principles density- 
functional theory. For free-standing superlattices, the 



ground state is the monoclinic Cm polar phase. Under 
the in-plane compressive epitaxial strain, it transforms 
to tetragonal polar Pimm phase, and under in-plane 
tensile strain, it transforms to orthorhombic Amml po- 
lar phase. All components of the static dielectric tensor 
(sij), the piezoelectric tensor (e;„), and the elastic com- 
pliance tensor (S^) were calculated as a function of the 
in-plane lattice parameter for 1/1 super lattice. The crit- 
ical behavior of some components of £,-j, e^, and 
tensors at the boundaries between different polar phases 
was observed. It was shown that the phase transitions 
between different polar phases are of the improper fer- 
roelastic type. The possibility of obtaining ultrahigh 
piezoelectric moduli using fine tuning of the epitaxial 
strain in superlattices was demonstrated. The compar- 
ison of the mixing enthalpy calculated for superlattices 
and disordered Bao.sSro.sTiOs solid solution modeled us- 
ing two special quasirandom structures SQS-4 revealed a 
tendency of the BaTi03~SrTi03 system to short-range 
ordering of cations and showed that short-period super- 
lattices are thermodynamically quite stable. 
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